function sse = FitDoubleGauss(p,xdata,ydata )
        FittedCurve = DoubleGauss(p,xdata);      
        ErrorVector = FittedCurve' - ydata;
        sse = sum(ErrorVector .^ 2);
end